First ground based measurement of atmospheric Cherenkov light from cosmic rays 
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A recently proposed novel technique for the detection of cosmic rays with arrays of Imaging 
Atmospheric Cherenkov Telescopes is applied to data from the High Energy Stereoscopic System 
(H.E.S.S.). The method relies on the ground based detection of Cherenkov light emitted from the 
primary particle prior to its first interaction in the atmosphere. The charge of the primary particle 
(Z) can be estimated from the intensity of this light, since it is proportional to Z^. Using H.E.S.S. 
data, an energy spectrum for cosmic-ray iron nuclei in the energy range 13-200 TeV is derived. 
The reconstructed spectrum is consistent with previous direct measurements and is one of the most 
precise so far in this energy range. 
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PACS numbers: 96.50.sb,98.70.Sa,96.50.sd 



I. INTRODUCTION 
A. Motivation 

Cosmic rays reach the earth at a rate of approximately 
1000 m~^. Their energy spectrum is steeply falling 
and remarkably featureless over ten orders of magnitude 
in energy The differential flux is well described by 
a power law ($ E^^) with a steepening of the spec- 
trum at a few PcV (the so-called "knee" ) and a flattening 
around 10 EeV (the so-called "ankle"). Despite advances 
in the field, the origin of cosmic rays is still unresolved. 
Supernova explosions are thought to be the major con- 
tributor at energies up to f PeV [1, 0], but conclusive 
proof is still missing. The High Energy Stereoscopic Sys- 
tem (H.E.S.S.) [1] has clearly identified supernova shock 
waves as sources of high-energy particles [5, 6]. However, 
the nature of these particles - electrons or cosmic-ray 
nucleons - remains under debate. 

The elemental composition of cosmic rays is similar to 
the composition of the solar system, if one accounts for 
propagation effects through the galaxy 0. At present 
the best measurements of elemental spectra in the en- 
ergy range 1 GeV to 0.5 PeV come from long duration 
balloon flights Q . Because of the decreasing flux of cos- 
mic rays and the limited collection area of these experi- 
ments («1 m^), it is hard to extend such measurements 
to higher energies. A further improvement in the ac- 
curacy and energy range of composition measurements 
could provide crucial information about the acceleration 
mechanism and propagation of these particles, and there- 
fore provide further clues about their origin. 

In 2001, Kieda et al. 9] proposed a new method for the 
measurement of cosmic rays with Imaging Atmospheric 
Cherenkov Telescopes (lACTs) . The central idea of this 
method is to detect the Cherenkov light emitted by pri- 
mary cosmic-ray particles (so-called Direct Cherenkov 
light) from the ground. While DC-light has been mea- 
sured in the past by balloon experiments [l^, |Tl[, the 
measurement from the ground takes advantage of the 
huge detection area («10^ m^) of lACTs, in principle 
enabling the extension of spectral and composition mea- 
surements up to ~ 1 PeV. Here we review this technique 
and describe its application to data from H.E.S.S. We 
present the measurement of the iron spectrum and give 
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FIG. 1: Schematic representation of the Cherenkov emission 
from a cosmic-ray primary particle and the light distribution 
on the ground and in the camera plane of an TACT. 

an outlook on future applications of this method. 



B. Technique 

When cosmic rays enter the atmosphere they emit 
Cherenkov light above an element-dependent energy 
threshold. The Cherenkov angle increaces with the den- 
sity of the surrounding medium. The emission angle of 
the DC-light therefore increases with increasing depth of 
the primary particle in the atmosphere, creating a light 
cone on the ground with a radius of roughly 100 m (see 
Fig. [1]). At a typical height of 30 km the particle interacts 
and a particle cascade is induced (Extensive Air Shower, 
EAS). The Cherenkov light from these secondary parti- 
cles creates a second, wider, light cone on the ground. 

The intensity of the DC-light is proportional to the 
square of the charge Z of the emitting particle, and can 
therefore be used to identify the primary particle. The 
challenge for detecting DC-light is to distinguish it from 
the much brighter EAS-light background (Fig. [2]). Be- 
cause the DC-light is emitted higher in the atmosphere, 
it is emitted at a smaller angle than the EAS-light, and 
is therefore imaged closer to the shower direction in the 
camera plane. A typical emission angle for DC-light is 
0.15° to 0.3°, whereas most of the EAS-light is emitted at 
angles greater 0.4° from the direction of the primary par- 
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FIG. 2: Simulated intensity distribution on the ground for 
the EAS-light and DC-hght of an individual 50 TeV iron nu- 
cleus, as a function of distance from the shower core, for two 
different first interaction heights (the shower core is defined 
as the intersection point of the shower axis on the ground). 
The zenith angle is 0°. The drop in DC-intensity at 100/120 
m reflects the first interaction height. The low intensity tail 
at larger radii is caused by Cherenkov light from fragments of 
the primary nucleus. 



tide (for a more detailed discussion see [9|). Cherenkov 
cameras, with pixel sizes of ~ 0.1° are therefore able 
to resolve the DC-emission as a single bright pixel be- 
tween the reconstructed shower direction and the center 
of gravity (cog) of the EAS-image in the camera plane 

(Fig. m- 

Unfortunately, the number of emitted DC-photons also 
depends on the emission height and on the energy of the 
primary particle (Fig. [3|). The height of first interaction 
of hadrons typically varies between 20 to 40 km, hence 
the total amount of emitted DC-photons of particles with 
the same atomic number and energy varies significantly. 
As will be discussed later, this makes the estimation of 
the primary charge more difficult. 

The energy range to which this technique can be ap- 
plied depends on the charge of the primary particle [9| . 
At lower energies the limiting factor is that the primary 
particle momentum must exceed the Cherenkov thresh- 
old. At very high energies, the EAS-light outshines the 
DC-light, making the detection of the latter impossible. 
The reason for this is that the intensity of the EAS-light 
increases approximately linearly with energy, whereas the 
amount of emitted DC-photons remains basically con- 
stant above a certain energy (see Fig. [31). 

The aim of the current work is to demonstrate that the 
technique of DC-light detection can be applied to instru- 
ments such as H.E.S.S. to measure the flux of cosmic-ray 
iron nuclei. Because of their large atomic number and 
high flux compared to other heavy elements, iron nuclei 
are well suited for DC-light detection. The lower energy 



FIG. 3: Total number of emitted DC-photons as a function 
of energy and first interaction height for an iron nucleus at 
a zenith angle of 0°, calculated using an atmospheric profile 
appropriate for the H.E.S.S. site. 



threshold for the detection of these nuclei is ^ 10 TeV. 

The H.E.S.S. instrument used for this measurement 
consists of four lACTs situated in the Khomas highland 
of Namibia, at a height of 1830 m above see level. Each 
telescope is equipped with a 960 pixel camera. Each pixel 
has an angular diameter of 0.16°, providing a total field 
of view of 5° diameter. The four telescopes are triggered 
in coincidence and image the Cherenkov light from EASs. 
The exact trigger conditions can be found in \1M ^^'^ 
are far below the applied analysis cuts described in the 
next section. As will also be described there, the proper- 
ties of the primary particle (such as direction or energy) 
can be reconstructed from the shower images. 

H.E.S.S. is a 7-ray experiment. The main challenge 
in detecting 7-rays is to distinguish them from the much 
larger background of hadronic cosmic rays (see for ex- 
ample [l3|). This background, recorded during normal 
7-ray observations, is now used to search for events with 
DC-hght. 



II. EVENT SELECTION AND 
RECONSTRUCTION 

A. Shower reconstruction and candidate event 
selection 

As for the standard H.E.S.S. analysis, the raw shower 
images are calibrated [14] and the pixel intensities are 
corrected to account for the loss of optical efhciency of the 
system over time [l^ . Afterwards the images are cleaned 
to remove low intensity substructure in the shower images 
and hence improve shower reconstruction. The image 
cleaning consists of a two-staged tail-cut, which requires 
pixels to have an intensity greater than 20 (10) photo- 
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electrons (pe) and a neighboring pixel with an intensity 
of 10 (20) pe. Afterwards, an "island cleaning" is ap- 
plied, where all pixels that are not connected to the pixel 
with the maximum intensity by neighboring pixels are 
removed. For a further reduction of shower fluctuations 
the same procedure is applied again with stronger cuts 
(200 and 100 pe). By default the soft-cleaned images are 
used for the shower reconstruction. The strong-cleaned 
images are only used instead if they contain more than 7 
pixels in more than two camera images. 

The showers are reconstructed using the standard 
stereoscopic Hillas analysis [l5| , whereby the shower (and 
therefore particle) direction and the intersection point of 
the shower axis on the ground are reconstructed by in- 
tersecting the major axes of the different shower images. 
The energy of the primary particles is reconstructed from 
the total image intensity /tot , the impact parameter i?corc 
(perpendicular distance from the shower axis to the tele- 
scope) and the zenith angle 0, by comparing these pa- 
rameters to simulations. The mean EAS-light yield at a 
fixed energy varies with the atomic number of the pri- 
mary particle, which introduces a systematic shift in the 
energy reconstruction between different elements. In this 
analysis all energies are reconstructed under the assump- 
tion that the primary particle is an iron nucleus. 

This energy reconstruction technique leads to a sys- 
tematic bias close to the energy threshold for detection 
[l3| . Therefore only events with a reconstructed energy 
greater than 13 TeV were considered, for which the en- 
ergy bias is less than 5% (the exact energy value is 12.59 
TeV, which corresponds to logig(i?/TeV) = 1.1 and will 
always be referred to as 13 TeV in the following). Addi- 
tionally, to avoid images truncated at the camera edge, 
only images with a center of gravity less than 2° from 
camera center are used. Finally, to select well recon- 
structed showers, only events that contain at least two 
camera images with an Aspect Ratio smaller than 0.75 
(Aspect Ratio = """^'^ width \ ^^^.^ considered. 

^ ^ image length, ' 



B. DC-light detection 

DC-Light can be identified as a single high intensity 
pixel between the reconstructed shower direction and the 
cog of the EAS-shower in the camera images (Fig. [ij. 
The main selection parameter for finding this DC-pixel 
in the camera image is the DC-ratio, defined as: 



Qdc 



I, 



max.ncighb. 
xcl 



(1) 



where /max.ncighb. is the maximum intensity of the neigh- 
boring pixels. The pixel with the minimum Qdc is de- 
termined in the relevant angular region in each camera 
image. The parameters used to constrain this region are 
the angular distance from the DC-pixel to the shower di- 
rection (AjJ,'^), to the cog of the EAS (Ao'g) and to the 
line connecting these two points (A^j;;). Afterwards a 
selection on the impact parameter is applied to ensure 
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TABLE I: Cut parameters for DC-pixel detection, 
nitions of the parameters see text. 
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FIG. 4: Mean DC-ratio Qdc as a function of the total inten- 
sity of the camera image Jtot for iron events from simulations 
(events which interacted in the atmosphere before passing the 
Cherenkov energy threshold were not considered since they 
contain no DC-light). The zenith angle is 0°. The error bars 
show the RMS of the distribution in each bin. The dotted 
line shows the selection function given in table IT] 



that the telescope is inside the DC-light cone. Finally, 
to avoid detector saturation effects, the intensity of the 
DC-light candidate pixel /dc-pixcI has to be below the 
saturation intensity. The exact cut values for the iden- 
tification of a DC-pixel were optimized using iron simu- 
lations and are summarized in Table IH To illustrate the 
applied cuts on one example, the distribution of Qdc for 
iron nuclei and the cut values are shown in Fig. 2] as 
a function of the total image intensity /tot- The mean 
Qdc value depends on /tot because at higher energies 
more and more Cherenkov light from secondary particles 
falls into the angular region of the DC-light, while the 
DC-intensity remains basically constant above a certain 
energy threshold. 

Once a pixel in the camera image fulfills all the men- 
tioned selection parameters the DC-light intensity /dc is 
reconstructed by subtracting the mean intensity of the 
neighboring pixels /ncighb. pixels from the DC-pixel inten- 
sity: 



-^DC — ^DC-pixel— < ^neighb. pixels > 



(2) 
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Z range 


representative element 


[10 * s ^sr ^ m ^] 


i-D 


P (Z=l) 


12. D 


7-9 


O (Z=8) 


1.43 


10-16 


Mg (Z=12) 


2.09 


17-24 


Ca (Z=20) 


0.56 


25-28 


Fe (Z=26) 


2.50 



TABLE II: Representative element and integrated flux be- 
tween 13 and 200 TeV, -Fi3°, of the reference composition 
Pi [HI for the five charge bands of the simulated fiux. 



III. DATA ANALYSIS AND SIMULATIONS 
A. Simulations 

As for all air shower experiments, the present anal- 
ysis relies on comparisons to Monte Carlo simulations. 
Simulations are used to calibrate the energy and charge 
fsection FlII Pp estimation and to determine the detection 
efficiency of the system (section|VB]) . In order to perform 
a direct comparison between simulations and data, sim- 
ulations have been produced for five different elements, 
representative of five different charge bands (see Table 
in|) . The charge bands cover a Z range from 1 to 28. 
The contribution of elements with Z > 28 is likely neg- 
ligible: although no flux measurements for these ultra 
heavy cosmic rays exist in the energy region of impor- 
tance here, measurements of energy spectra for elements 
with Z > 28 at lower energies P, [l3| and flux predictions 
in the TeV energy region are more than three orders 
of magnitude below the flux of the iron band. 

The contribution of each charge band to the total simu- 
lated flux and the energy distribution of the events inside 
each band are weighted according to the measured fluxes 
given in and [H (referred to in the following as ref- 
erence composition). Ref. [l^ shows a parameterization 
of the different elemental spectra, obtained by combin- 
ing measurements from several experiments. The errors 
given on the absolute flux normalization are <10%. How- 
ever, they probably do not reflect the entire uncertainty 
in the spectra, since systematic uncertainties may have 
been underestimated in the individual data sets [16| . Ad- 
ditionally, in the energy range of interest here, specifically 
for the calcium band, the values given are extrapolated 
from measurements at lower energies. Ref. [l| presents 
the same parameterization with more recent data for the 
proton, helium and iron flux. The normalization of the 
flux for these elements differs by approximately 25% be- 
tween the two parameterizations. For the mentioned rea- 
sons we assume this difference to be a realistic error in the 
integral fluxes of the different charge bands between 13 
and 200 TeV. In the following comparisons between the 
data and the simulations, this error is always included. 
Since the 25% error is still somewhat arbitrary, we will 
also discuss the effect of a more conservative error of 50% 
on the presented measurements at the end of this work. 

The shower parameters of relevance here depend on 
the details of high energy hadronic interactions. To as- 



sess the systematic errors arising from uncertainties in 
these interactions, the analysis is performed with sim- 
ulations based on two independent hadronic interaction 
models, SIBYLL 2.1 [1^ and QGSJET Olf ^20^. (The 
newer version QGSJET 02 [2lj was not available at the 
time of the analysis. However, the hadronic interaction 
uncertainty estimation with QGSJET Olf should be more 
conservative, since more recent hadronic models are ex- 
pected to model the interactions more accurately). The 
simulations for both models were performed using the 
shower simulation program CORSIKA 6.0321 22]. For 
each model a total of ~ 10^ showers were simulated in an 
energy range from 1 to 200 TeV. The zenith angle of the 
simulations was chosen to match the mean zenith angle 
of the data set (13.6°) described in the next section. 

B. Data 

The data considered here were taken between 2004 and 
2006 with the full four telescope H.E.S.S. array. Because 
cosmic rays arc deflected by magnetic fields in the galaxy, 
their flux in the measured energy region is expected to 
be very close to isotropic |7|- Therefore it is possible to 
use all available H.E.S.S data, independent of the target 
position. However, to reduce systematic uncertainties 
due to zenith angle effects on detection efficiencies and 
energy reconstruction, only data runs with a mean zenith 
angle smaller than 22° are considered. After standard 
quality selection criteria and dead time correction, the 
data set amounts to 357 hours of observation time. 

In total, 35364 events in the energy region from 13 to 
200 TeV passed the selection criteria. One example of an 
event with DC- light in all four telescopes is shown in Fig. 
[S] High intensity pixels close to the reconstructed shower 
direction are evident in all four images (indicated by ar- 
rows). In agreement with the expected trend from Fig. 
[21 the highest /lowest DC-pixel intensity also corresponds 
to the largest/smallest impact parameter, respectively. 

Fig. [H] shows the DC-light detection rate as a function 
of the telescope multiplicity iVtci (number of telescopes in 
which DC-light is detected simultaneously) for data and 
simulations for both hadronic models. While the shape 
of the distributions agrees well, the event rate is higher 
by «25% for the SIBYLL simulations. This difference 
between the models gives an estimate of the systematic 
error introduced in the analysis due to hadronic inter- 
action uncertainties. Taking this systematic error into 
account, the simulated rates and the data are consistent. 



C. Background 

The detection rate of events with DC-light is ex- 
pected to have some background due to misidentifica- 
tions. These misidentifications can occur due to shower 
fluctuations, which can lead to single high intensity pixels 
in the EAS-light images. The rate of false detections can 
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FIG. 5: A measured event with indications of DC-light in all four cameras images (indicated by arrows), after high threshold 
image cleaning. The reconstructed shower direction is shown by a cross (x) in each image. The reconstructed energy of this 
event is 50/48 TeV based on QGSJET/SIBYLL simulations. The reconstructed impact parameter and DC-hght intensity for 
each telescope are shown in the lower panels in each image. The energy and impact parameter resolutions are «20% and «20 
m, respectively. The white points mark disabled pixels. 



be estimated using proton simulations. Protons domi- 
nate the cosmic-ray flux in the energy region of interest 
(see Table lll| . Additionally, protons emit only a negligi- 
ble amount of DC-light compared to their EAS-light yield 
at these energies [9,], therefore any detection of DC-light 
in proton simulations can be considered a fluctuation of 
the EAS-light yield and therefore a false detection. As 
shown in Fig. [6l the expected misidentification rate for 
the proton band is «10% of the measured event rate in 
the data for A^tci^l- However, for iVtci >2 there are no 
misidentified events found in the proton band for either 
hadronic interaction model. The upper limit derived for 
the misidentification rate for these events is almost two 
orders of magnitude below the detected event rate in the 
data. Events with iVtoi >2 can therefore be considered 
as essentially background-free. 



background estimation in the presented analysis, only 
events with A'tci >2 were considered. They will be ref- 
ered to in the following as DC-events. In total, 1899 
DC-events were found in the data. The resolution of 
the shower parameter reconstruction for these events is 
w0.1° for the shower direction, Ki20 m on the shower core 
position and ^15% on the primary energy. 



D. Primary charge reconstruction 

The elemental composition of the DC-events can be 
estimated using the Z dependence of the DC- light inten- 
sity. The reconstructed charge Z* is defined as: 



In order to minimize systematic uncertainties due to 



z* = d{E,e)^iDc 



(3) 
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FIG. 6: Event rate for different telescope multiplicities A^toi 
compared to the expected rate from the simulation of the 
cosmic-ray flux using the SIBYLL and QGSJET interac- 
tion models. Additionally, the background event rate from 
the proton band, estimated using SIBYLL simulations, is 
shown. No background events were identified for A'td > 2 
and upper limits are shown (this result is the same for both 
hadronic models, which is why only the SIBYLL simulations 
are shown). An upper limit at 95% confidence level is shown 
for A'tei = 2. The upper limits for A'td > 2 are derived under 
the assumption that they have the same ratio to the detected 
event rate as the Atei = 2 events. This is a conservative es- 
timate since the fraction of misidentification is expected to 
decrease with increasing Ntai- 



where d{E, 9) is a factor that normalizes the mean of 
the Z* distribution from iron simulations to the atomic 
number Z of iron. The energy dependence of d is due 
to the energy dependence on the number of emitted DC- 
photons, since the emission for iron nuclei is not satu- 
rated in the lower part of the observed energy range (Fig. 
[3]). The zenith angle dependence of d arises from the 
increasing distance between the average first interaction 
point to the telescopes with increasing 6. 

The charge resolution achieved using Z* is energy de- 
pendent and improves for higher energies, for two prin- 
cipal reasons: 

1. The separation of the DC-light intensity distribu- 
tions, and hence the Z* distributions, for differ- 
ent elements is maximized when the elements com- 
pared have high enough energies that their DC- 
light emission is saturated. The saturation energy 
increases with charge. The heaviest element in this 
analysis is iron, for which the saturation energy is 
«50 TeV, which means that the charge separation 
of Z* continues to improve up to this energy. 

2. A significant fraction of the detected DC-events 
are dominated by emission from secondary particles 
created in the first interaction. The reconstructed 
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FIG. 7: Mean reconstructed charge Z* as a function of the 
true charge Z for DG-events in an energy range of 1.5 < 
logjg(i?/TeV) < 1.7 for both hadronic models. For clarity 
the X-axis is shifted by -1-1 for SIBYLL. The error bars show 
the RMS of the distribution in each bin. 



charge Z* for these events is lower than the charge 
of the primary particle. Simulations show that the 
fraction of DC-events dominated by the DC-light 
from secondary particles decreases with energy. For 
iron it is «60% at low energies (wlS TeV) and 
drops to «35% above the saturation energy. 

Fig. [7] shows the charge resolution obtained in an in- 
termediate energy range for both hadronic models. The 
charge resolution achieved is not sufficient to assign the 
charge of the primary on an event-wise basis to one of 
the four charge bands. However, as will be shown later, 
it is possible to measure the fraction of elements belong- 
ing to the iron band in the data on a statistical basis, 
and therefore estimate the iron flux. 

The main reason for the relatively broad distribution 
of Z* for each element is that the DC-light intensity de- 
pends not only on the charge of the primary, but also 
on the emission height of the DC-light (Fig. [3]). The 
mean emission height is determined by the first interac- 
tion height distribution of the primary particles in the 
atmosphere. As mentioned this varies significantly from 
event to event, with a FWHM of «10 km. This depen- 
dence on the emission height is one of the reasons for 
the observed bias in the charge reconstruction at lower 
atomic numbers, since the mean first interaction height 
varies between elements. A second reason is that the re- 
construction of Z* is normalized to the charge of iron, 
which has not saturated its DC-emission in this energy 
region. 
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In principle a more accurate measurement of the pri- 
mary charge could be achieved by constraining the mean 
emission height ft-DC- This is possible because hjnc is 
directly related to and i?corc via: 

hue - (4) 

^DC 

This equation follows directly from the geometry of the 
emission, where A^^ is the mean Cherenkov angle un- 
der which the DC-light is emitted. However, in the case 
of the H.E.S.S. telescope system, little is gained by in- 
cluding the /iDC dependence in Z*, since the pixel size 
of 0.16° and the spread of the direction and impact pa- 
rameter reconstruction for DC-events are too coarse to 
provide a sufficiently precise determination of A^^ and 
_Rcorc- It should be noted that the limited reconstruction 
accuracy of i?core also limits other techniques, such as 
that proposed in which rely on this parameter. 



The larger difference between the two hadronic models 
for Xinax compared to the other parameters shown is not 
completely unexpected. The shower maximum depends 
on the exact modeling of the fragmentation processes in 
the EAS-shower. The differences in the fragmentation 
process between hadronic models are well known (see for 
example [l^, [ISl)- In contrast to the shower maximum, 
the distributions of i?corc, ^dc of the DC-events 

are dominated by the properties of the DC-light. These 
are easier to model because they are completely deter- 
mined by the distribution of the first interaction height 
for a given energy and atmospheric profile. Since com- 
position measurements with the DC-light technique rely 
primarily on the DC-light intensity, they are expected to 
be relatively model independent. 



V. SPECTRUM EXTRACTION 
A. Iron fraction 



IV. SYSTEMATIC CHECKS 

Since the reconstruction of the energy spectrum of 
cosmic-ray iron relies on Monte Carlo simulations, it is 
important to demonstrate adequate agreement between 
the measured and simulated distributions of parameters 
used in the reconstruction. Fig. [8] shows such compar- 
isons for the impact parameter, the DC-light intensity, 
the mean Cherenkov angle and the shower maximum 
^max- The shower maximum is the atmospheric depth 
at which the maximum number of Cherenkov photons is 
emitted in the EAS. The panels beneath the distribu- 
tions show the mean values of the distributions for the 
data and both hadronic models. The error bars on the 
simulated points include both the statistical uncertainties 
and the uncertainty in the cosmic-ray mass composition 
in this energy regime. 

The distributions of i?corc, ^dc and A^'J, show a good 
agreement between data and simulations. Their mean 
values agree within 1 a. For the height of the shower 
maximum Xmax, a shift of ~5% between SIBYLL simu- 
lations and the data is apparent. However, no significant 
shift in this parameter is present in the QGSJET simu- 
lations. This difference between the models is again an 
estimate for the systematic error arising from hadronic 
interaction uncertainties. Within this systematic uncer- 
tainty data and simulations agree reasonably well for 
^max- However, that the mean Xmax values for both 
hadronic models are larger than the data is an indication 
that there might be additional systematic uncertainties 
in the reconstruction of this parameter. These could for 
example come from the uncertainty in the atmospheric 
profile at the H.E.S.S. site, which is estimated to be ~3 
g cm~^. However, since no cut is applied on X^ax, these 
systematic effect could only affect the results indirectly 
and must be smaller than the hadronic model uncertain- 
ties. 



The first step in the derivation of the flux of iron nu- 
clei is the measurement of fraction fcpc of iron events 
among the DC-events, k-pc is estimated via a fit of a 
two-component model to the Z* distribution of the data. 
The first component of this model is the Z* distribution 
of simulated iron nuclei. The second component is a sum 
of the Z* distribution of lighter nuclei. The relative com- 
position of the lighter charge bands (= all except the iron 
band) is kept fixed to the reference composition, so that 
/cpc is the single free parameter of the fit. 

The fit was performed in five energy bands. The Z* 
distribution from the data and the fitted model using 
SIBYLL simulations are shown in Fig. [9l The fit results 
(/cpc) and corresponding values (Xdata) i^or the fits are 
summarized in the Table [m] for QGSJET and SIBYLL. 
The values for fcpc for the two hadronic models agree with 
each other within statistical errors Afcgt (the standard 
deviation of the fit result) for individual energy bands. 
However, the QGSJET values are shifted by ps-O.I. 

Apart from statistical uncertainties, fcpc is affected by 
the systematic uncertainty in the assumed composition 
of the lighter nuclei Afccomp- This is estimated by varying 
the weight of the individual lighter charge bands of the 
fitted model by 25% and performing the fit for each possi- 
ble combination. The minimum and maximum deviation 
between these fit results and the previously obtained /cpo 
value are then taken as errors on the composition uncer- 
tainty. Since these errors are close to symmetric, their 
absolute values are averaged to give Afccomp- 

Table IIIII also shows the expected iron fractions kp^ 
from simulations of the cosmic-ray flux assuming the ref- 
erence composition. These values agree with the best 
fit values. This implies, together with the reasonable 

/ndf values of the fits, that no significant deviation 
from the reference composition can be found in the data. 
For the iron band this statement will be quantified in the 
next section. 
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FIG. 8: Comparison of shower parameter distributions for data and simulations of the cosmic-ray flux. The bigger panels show 
normalized distributions of Rcore, /dc, Adc and Xmax for SIBYLL simulations and data. The smaller panels underneath the 

distributions show their mean values and the mean value for the different charge bands of which the simulated flux is composed. 
Additionally, they show the mean value for the distributions for the QGSJET simulations. When comparing the mean values 
for the simulated fluxes one should bear in mind that roughly equal contributions to the error bars come from statistical errors 
and from uncertainties in the reference composition, which is the same for both models. The systematic difference between the 
mean values for the shower maximum Xmax is not unexpected, since this quantity is difficult to treat in the simulations (see 
text). 
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logio(E/TeV) 


fcpg i AA)fit i ^^comp 


Xdata /ndf 


I.ref 




1.1-1.3 


0.56±0.047±0.026 


6.9/15 


0.55±0.09 




1.3-1.5 


0.64±0.049±0.027 


31.5/17 


0.64±0.08 


SIBY 


1.5-1.7 


0.77±0.054±0.019 


15.4/15 


0.70±0.07 


1.7-1.9 


0.66±0.097±0.013 


12.0/14 


0.80±0.06 




1.9-2.3 


0.93-h0.07-0.151±0.008 


3.4/12 


0.84±0.04 


H 


1.1-1.3 


0.47±0.050±0.026 


11.6/15 


0.50±0.10 


H 


1.3-1.5 


0.55±0.059±0.029 


13.3/17 


0.66±0.08 


1-5 

tn 


1.5-1.7 


0.70±0.063±0.023 


14.3/15 


0.75±0.06 


O 
C? 


1.7-1.9 


0.54±0.128±0.034 


9.3/14 


0.79±0.06 


1.9-2.3 


0.69±0.160±0.012 


6.2/12 


0.84±0.05 



TABLE III: The best fit value of the iron fraction in the data fepe and the values of the fit are shown for both hadronic 
models in five energy bands. The error of fcpe is composed of the statistical error of the fit Afcet and the systematic error from 
the uncertainty in the assumed composition of the lighter nuclei Afccomp (see text). Additionally shown for comparison is the 
fraction of iron in the simulated cosmic-ray flux (kpg ) for both models. 
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FIG. 9: Z* distributions of data (open circles) and the fit- 
ted SIBYLL model (black boxes) in five energy bands. The 
curves show the charge distributions of iron (dotted) and of 
the lighter charge bands (dashed) that compose the fitted 
model. The fit results are summarized in Table lllTl Note that, 
as mentioned in section IlII Dl the charge resolution improves 
with energy, observable by a narrowing of the Z* distribution 
of iron with increasing energy. 



B. Iron Flux 

Since the identification of DC-events is effectively 
background-free, the differential iron fiux 4i{E) can be 
estimated as: 



HE) 



Nuc{E) 



A,b{E) -AE-t 



■ k 



Fg, 



(5) 



where Nj^ciE) is the number of detected DC-events in 
the energy interval from E to E + AE, t is the total live- 
time of the data-set and Acs is the mean effective area 
times the field of view of the detector, averaged over the 
zenith angle of the observations, taking into account the 
efficiency of selection cuts. Acs is derived from simula- 
tions of iron nuclei via: 



Ar 



jV^C(^).^MC.^MC 



(6) 



where N^'~^{E) is the total number of simulated events 
in the energy interval from E to E + AE and is 
the corresponding number of identified DC-events. 
and 17^^ are the area and angular region over which the 
simulations were performed. 

The energy spectrum is measured in the five energy 
bins of the fcpo fit. The result is shown in Fig. [10] for 
both hadronic models together with the highest energy 
baloon measurements. The derived spectrum agrees well 
with these measurements for both models. The measured 
spectrum is fitted well by a power law 4>{E) — 4'o{^^y)~'' ■ 
The best fit values for the SIBYLL spectrum are given by 
(/)o = (0.029 ±0.011) m-^sr-i TeV^i and 7 = 2.76 ±0.11 
with an xV^df of 3.0/3. For the QGSJET spectrum the 
best fit values are 0o = (0.022 ± 0.009) m^^gj.-! TeV"^ 
and 7 = 2.62±0.11 with xVndf of 5.3/3. The integrated 
fiux above 13 TeV is F{> 13TeV) = (1.9 ± 0.7) • lO"'' 
s-^sr-i m-i for SIBYLL and F(> 13TeV) = (2.3 ±0.9) • 
10--* s-^sr-^m-i for QGSJET. 

Since both spectra are derived using the same data- 
set, the differences in the spectral index A7 = 0.14 and 
integrated flux AF/F = 17% again provide an estimate 
of the systematic error due to hadronic interaction un- 
certainties. Additional systematic errors, arising from 
uncertainties in the atmospheric profile and the absolute 
detection efficiency of the H.E.S.S. instrument are dis- 
cussed in detail in [l^ and lead to a systematic error of 
20% in the integrated flux and A7 = 0.1 in the spectral 
index. The effect of the systematic error Afccomp in fcpe on 
the spectrum amounts to A7 = 0.015 and AF/F = 5%. 
Assuming a more conservative error of 50% in the inte- 
gral fluxes of the lighter elements in the model of the fepc 
fit increases A/ccomp by 0.03 on average and leads to er- 
rors of A7 = 0.04 and AF/F = 11% in the presented 
spectrum. This error is still small compared to the pre- 
viously mentioned uncertainties. The total systematic 
uncertainty of the measurement is therefore estimated to 
A7 = 0.17 for the spectral index and AF/F = 28% for 
the integrated flux. 
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FIG. 10: DifFerential iron energy spectrum measured with H.E.S.S. for the hadronic models QGSJET and SIBYLL multiplied 
by E^'^ for better visibility of structures. The spectral points for both models are measured for the same energies. For better 
visibility the SIBYLL points were shifted 10% upwards in energy. The error bars show the statistical errors. The systematic 
flux error in each bin is 20%. The measurements from balloon experiments with data points at the highest energies are shown 
for comparison [25l , [2^ (a compilation with more measurements from balloon experiments and space born measurements can 
be found in [l[). For a better visibility no horizontal bars marking the bin ranges are shown, they can be found in the respective 
papers. When comparing the measurements one should bear in mind that the experiments have different charge thresholds for 
their definition of the iron band (see legend). 



The statistical error on the measured iron flux is com- 
parable to these systematic errors. This means that with- 
out an improvement in the latter, the total error of the 
measurement can not be significantly reduced by increas- 
ing the exposure time of the data set. However, an in- 
creased data-set would enable one to extend the measure- 
ment towards higher energies. We note that, despite the 
systematic uncertainties, the iron flux determined with 
this technique is one of the best measurements in this 
energy range. The good agreement between the mea- 
sured fluxes from balloon experiments and those given 
here lends confidence to the results from both techniques. 



VI. SUMMARY & OUTLOOK 

A technique for the detection of cosmic rays by resolv- 
ing the Cherenkov emission from primary particles has 
been presented and applied to H.E.S.S. data. As a result 
1899 events with Direct Cherenkov light in at least two 
telescopes were detected and it was shown that these DC- 
events can be considered as background-free. Different 
parameter distributions of these events were compared 
to simulations using two different hadronic interaction 
models and good agreement with the data was found for 
both. The strong correlation between the DC-light and 
the charge of the primary shower particle made a charge 
estimate possible, from which the energy dependent frac- 
tion of iron in the data was derived. The energy spectrum 
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of iron nuclei was determined in an energy range of 13 to 
200 TeV. The result confirms the flux measurements from 
balloon experiments with an independent technique and 
is one of the most precise measurements in this energy 
range. 

Future improvements of the DC-light technique could 
extend the energy range of the measurement to an energy 
of ~ 1 PeV. Besides larger statistics, this extension re- 
quires additional separation power of the DC-light from 
the EAS-light. The reason for this is that the DC-hght 
yield remains constant above a certain energy while the 
EAS light yield increases approximately linearly with en- 
ergy. As shown in , additional separation power can be 
achieved using the time structure of the DC-light, since 
it arrives with a typical delay of 4 ns with respect to 
the EAS light. This fact could not be exploited in the 
analysis presented because the H.E.S.S. data used here 
were taken with the standard integration window of 16 
ns. However, current and planned Cherenkov telescopes, 
which routinely store pulse timing information [13, , 
may take advantage of this characteristic. 

Due to the strong dependence of the DC-light yield 
on the charge of the primary particle, the DC-light tech- 
nique has great potential for composition measurements. 
The limiting factor is currently the accuracy of shower 
reconstruction, needed to constrain the emission height 
of the DC- light. Since the typical shower images in the 
present work contain ^100 pixels, the limitation in the 
reconstruction accuracy of the shower with the Hillas 
technique arises from the strong fluctuations in hadronic 
showers and not from the limited angular resolution of 
the system. This puts a physical limit to the charge res- 
olution when the showers are reconstructed using this 
technique. 

In order to reduce the sensitivity to shower fluctua- 
tions, the DC-light itself could in principle be used to 
reconstruct the shower. One simple extension of the cur- 
rent method would be to determine the principal axes of 
the shower images as the line connecting the center of 
gravity of the image to the DC-pixel. However, exploita- 
tion of this technique would require pixels with smaller 



angular scale for the accurate localization of the DC-light 
spot. Another possibility, with an array of many nearby 
telescopes, would be to reconstruct the shower from the 
DC-light intensity distribution on the ground. Both tech- 
niques would take advantage of the very small fluctua- 
tions of the DC-emission. A quantitative statement on 
the level of improvement for both techniques would re- 
quire detailed simulations and is beyond the scope of this 
paper. 

Finally, the agreement of the distribution of the recon- 
structed charge Z* for large charges [Z* > 28) disfavors 
a significant contribution from ultra heavy elements to 
the cosmic-ray flux, as expected. It should also be noted 
that the present analysis is not sensitive to exotic states 
of matter, such as quark matter or magnetic monopoles 
as proposed in [l^. The high charges of these states 
{Z ^ 100) imply DC-intensities which would saturate 
the photomultipliers in the H.E.S.S. cameras. 
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